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ABSTRACT 

We used data from the Helioseismic and Magnetic Imager (HMI), and Atmospheric Imaging As- 
sembly (AIA) on the Solar Dynamics Observatory (SDO) to study coronal loops at small scales, 
emerging in the quiet Sun. With HMI line-of-sight magnetograms, we derive the integrated and un- 
signed photospheric magnetic flux at the loop footpoints in the photosphere. These loops are bright 
in the EUV channels of AIA. Using the six AIA EUV filters, we construct the differential emission 
measure (DEM) in the temperature range 5.7 — 6.5 in log T (K) for several hours of observations. 
The observed DEMs have a peak distribution around log T ~ 6.3, falling rapidly at higher temper- 
atures. For log T < 6.3, DEMs are comparable to their peak values within an order of magnitude. 
The emission weighted temperature is calculated, and its time variations are compared with those 
of magnetic flux. We present two possibilities for explaining the observed DEMs and temperatures 
variations, (a) Assuming the observed loops are comprised of hundred thin strands with certain radius 
and length, we tested three time-dependent heating models and compared the resulting DEMs and 
temperatures with the observed quantities. This modeling used Enthalpy-based Thermal Evolution 
of Loops (EBTEL), a zero-dimensional (OD) hydrodynamic code. The comparisons suggest that a 
medium frequency heating model with a population of different heating amplitudes can roughly re- 
produce the observations, (b) We also consider a loop model with steady heating and non-uniform 
cross-section of the loop along its length, and find that this model can also reproduce the observed 
DEMs, provided the loop expansion factor 7^5- 10. More observational constraints are required to 
better understand the nature of coronal heating in the short emerging loops on the quiet Sun. 
Subject headings: Sun: photosphere — Sun: surface magnetism — Sun: atmosphere — Sun: corona 



1. INTRODUCTION 

A part of the magnetic field originating in the pho- 
tospheric sub-surface layers reaches higher up in the 
solar atmosphere and forms loop like structures, the 
building blocks of solar corona. These loops harbor 
plasma, which is heated up to a few million Kelvin, 
much higher than the photospheric temperature. Find- 
ing the source and nature of energy required to heat 
the corona along with the process of heating is one of 
the most sought after questions in the field of astro- 
ph ysics (for reviews on coronal heating, see for exam- 
ple lZirkerlll993l: iNarain fc Ulmschneiderl 119961 : iKlimchukl 
120061 : lRealell2010n . Studying the dynamics of the plasma 
filled loops is important to understand the heating mech- 
anisms responsible for these high temperatures. Obser- 
vational, theoretical, and numerical advances have been 
made over several decades to understand the physics in- 
volved in these processes. Some of the early works on 
this subject include the ideas of damping of magnetohy- 
drodynamic waves in the lower corona to heat the solar 



a stea dy- st ate heating process. 



Parker! 



atmosphere (lAlfven ) [l947l for a recent review on waves 
i solar corona, see INakaria kov fc Verwichte1l2005f ). 
From the early x-ray observations ( Vaiana et al.lll973f ). 



it became evident that the solar corona is confined in the 
form of loops outline d by the underlying photospheric 
magnetic field. Later. iRosner et all (|1978t l. gave an an- 
alytical model for the quiescent coronal loops, assum- 
ing that these structures are in hydrostatic equilibrium. 
They suggested that the observations are indicative of 
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(|1994f h and iCargill fc Klimchuld (|1997h put forward the 
idea of intermittent and impulsive (nanoflare) heating, 
as a viable mechanism. It is now generally believed, and 
widely accepted that the magnetic field plays an impor- 
tant role in generating and transporting the energy re- 
quired to maintain the temperatures of the corona. It 
remains unclear and difficult to identify the dominant 
process responsible for heating of the solar atmosphere. 

As the diagnostics of tenuous coronal plasma improved 
with the advent of high spatial and temporal resolution 
space based instruments, an alternate but relevant de- 
bate emerged within the community, namely, the fre- 
quency of required heating events. The plasma filled in 
the loops respond to the impulse of heating, and this 
depends on whether the plasma is reheated before it is 
completely cooled down (high frequency model — steady 
heating), or not (low frequency model — nanoflares). 
Should either of these models operate, they predict cer- 
tain physical properties of the lo ops, which can be com- 
pared with the observations (sec Reale 2010, for a broad 
review on coronal loop observations and modeling). 

With a wide range of field strengths and sizes of mag- 
netic elements, coronal loops also have wide tempera- 
ture and length distributions. Usually the loops are 
classified as "hot" (T > 2 - 3 MK), and "warm"(T w 
1 — 2 MK) depending on their temperature regime. 
Both steady and impulsive heating models have been 
extensively used to explain the observed temperatures, 
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loop intensity structure e.t.c. Studies indicate that 
the hot pl asma is consistent with both steady heat - 
ing models ([Warren et alj|2010t iWinebarger et al.ll2011l) 
and impulsive heating models (|Tripathi et al.l [20101 : 
IViall fc Klimchul [2012]) . The warm loops are found 
to be continu ously evolving and not i n equilibrium (see 
for example, lUgarte-Urra et all I2009D , and their prop- 
erti es are well explaine d by impulsive heating mod- 
els ([Spadaro et al.ll2003() . It is also suggested that the 
age of an active region might play an important role 
in det ermining the dominance of one process over the 
other ([Schmelz fe Pathakl l20l"l lUgarte-Urra fe Warrenl 
[20121) . 

The active regions are well studied both in terms of ob- 
servations and modeling. However, the situation is not so 
clear in the case of small loops in the quiet Sun. The clas- 
sification of "hot" and "warm" loops may not be relevant 
in these features, owing to their compact magnetic struc- 
ture and narrow temperature range compared to the ac- 
tive regions. These short loops are connected to magnetic 
bipoles in the photosphere. Their origin can be traced to 
either flux emergence, or convergence of opposite polar- 
ities with reconnection. The magnetic fluxes associated 
with these regions are typically in the range of 10 19 — 10 20 
Mx. The electron number density in such loops, mea- 
sured using density sensitive lines is in the order of 10 9 
cm" 3 (lUgarte-Urra et aIll2005t[Plrez-Suarez et~aTll2008t 
iDoschek et al.ll2010D . 

In this study, we are primarily interested in under- 
standing the nature of the heating that produces the ob- 
served 1 — 2 MK temperature in these small bipoles, in 
particular the frequency of heating events. Also, to bet- 
ter understand the relation between photospheric mag- 
netic field and the coronal loop temperatures, we chose 
to study emerging flux events. In these events it is easy 
to identify the loops, and their footpoints in the pho- 
tosphere. We follow their formation and evolution over 
many hours. In the following section, we present the 
observational results. Section [5] describes the loop mod- 
eling and the simple heating models we tested in this 
work. Finally, we summarize the results, and discuss 
some relevant aspects that require further investigation. 

2. OBSERVATIONAL RESULTS 

In this section we give a brief note on the datasets used, 
and present the results derived, namely, the photospheric 
magnetic flux, and coronal temperatures. The line-of- 
sight magnetograms ob served with the Helioseismic and 
Magn etic Imager (HMI. lScherrer et al.ll2012tlSchou et al.l 
|2012|) . and the intensity images from the E UV channels 
of At mospheric Imaging Assembly (AIA, iLemen et al.l 
are used. HMI and AIA are two of the three in- 
struments onboard Solar Dynamics Observatory (SDO, 
Pesncll ct al. 2012). Data are taken from 2011 February 
10, and 2012 March 17 observations, spanning for about 
12 hr each. SDO observes the full disk of the Sun con- 
tinuously in different filters with a high cadence of 12 s. 
We selected a region near disk center with the criteria 
that, we see emerging magnetic field and coronal loops 
close to the beginning of the selected time sequence. A 
few cases of evolved bipoles are also considered. 

AIA data contain time sequences from 94 A, 131 A, 
171 A, 193 A, 211 A, and 335 A EUV channels. Data 
are processed with standard procedures available in the 
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Fig. 1. — A context image showing one of the analyzed bipoles 
as seen in AIA 193 A channel. The regions marked with black 
(near loop top), and white (background regionjboxes are used for 
further analysis to produce DEMs (see Figure __a).) 

solarsoft library. Alignment between the data from all 
these channels is crucial. Using 171 A images as refer- 
ence, and cross-correlation technique, we aligned all data 
to within a pixel. The emerging bipoles are identified 
both in HMI, and AIA. The tracked data cubes of such 
bipoles are extracted for further analysis. To enhance the 
signal-to-noise ratio, we prepare the 12 s cadence AIA 
data to 1 minute cadence by averaging five exposures in 
each channel. Next, to derive the physical properties of 
the plasma, we adopt the differential emission measure 
(DEM), which is related to the electron number density 
(n e ), and the line-of-sight plasma temperature gradient, 
and defined as 

, dh 



dT' 
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We use data from six AIA EUV channels, 
along with the filter response^ as input to con- 
struct DEM(T) (cm" 5 K" 1 ), at each pix el us- 
ing xrt_d.em_ite rative2.pro ([Golub et al.l 120041 : 
Web er et al.l 12004 distributed in solarsoft). In this 
program, initial DEM is guessed and folded through the 
filter responses to generate model observations, which 
are iteratively used to reduce the % 2 between the original 
and modeled observations. Thi s program uses a much 
tested IDL routine mpfit.pro (Markwardt 2009), that 
performs a Levenberg-Marquardt technique to solve the 
least-squares problem. 

We show results from four emerging bipoles in this 
work. In Figure __ we plot the region-of-interest for one 
of the bipoles analyzed, as seen in AIA 193 A. The im- 
age saturates at 750 DN s _1 . For these bipoles, a series 
of DEMs are constructed near the loop top (for example 
from a region with in the black box shown in Figure [[J) at 
each pixel in a 6" x 6" region, over several hours of obser- 
vations. The predicted intensities from forward modeling 
of the derived DEMs match the observed intensities with 
in the limits of errors. Since we restrict the DEMs within 
a limited range of temperature, the predicted intensities 
will be lower limits of the observed values. In Figure __ 
we plot the average emission from this area as a function 
of temperature (log T), for all times. Each panel cor- 
responds to a bipole. The dots denote the time depen- 
dence. At any given temperature, to show the emission 

The filter responses of 94 and 131 A channels are empiri- 
cally modified to include contributions from Fe IX and Fe XII for 
94 A, and from Fe VIII and Fe XI for 131 A. The revised response 
functions can be obtained using aia_get_response with a keyword 
chiantif ix, available in solarsoft. 
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Fig. 2. — DEMs of four bipoles obtained from the observations 
plotted as a function of temperature for all times (dots). A small 
offset in temperature is given to the DEMs to show the temporal 
distribution. Solid line is the temporal median of DEMs obtained 
for the respective cases. Observed DEMs have a peak around log T 
of 6.2 — 6.3. Dashed histogram in panel (a) is the temporal median 
of five DEMs, at random times, obtained from a background region 
close to the corresponding bipolc. 

distribution in time, we gave a small offset to DEMs in 
temperature (and that is the reason we see a small spread 
of DEMs along log T) . Additionally, the temporal distri- 
butions also give a sense for the errors in the DEMs. The 
solid lines are respective temporal medians for all DEMs. 
They have a peak close to log T (K) of 6.2 — 6.3. At 
higher T, they show a rapid decline and also the DEMs 
are not well constrained. On the other hand, at lower 
T, the emission stays comparable to the peak emission. 
Similar results were obtained using Hinode /EIS observa- 
tions, but for a coro nal hole bright point (c.f. Fig. 12, 
iDoschek et aIll20Toh . 

We also note that there is no background subtraction 
to the data in DEM analysis. The small loop structures 
we analyzed have their loop apex, and footpoints in the 
same plane along the line-of-sight, much of the emission 
contribution may be primarily dominated by the loop 
apex with a part of it originating from the footpoints. 
To compare the background contribution to the resulting 
DEMs, for example (a), we considered a 10 x 10 pixel 
background region adjacent to that bipole (marked with 
a white box in Figure [T]). The DEMs are constructed 
for this region at five random times, and the temporal 
median is plotted as a dashed histogram in Figure EJa). 
This shows that the observed DEM lies well above the 
background, not only for log T in the range 6.2 — 6.4, but 
also at low temperatures (log T < 5.8). The small scatter 
in the DEM at low T suggest that this emission is real 
and is not an artifact of the DEM inversion. The DEMs 
are now used to derive an emission weighted temperature 
{Tdem) using the relation 

S {DEM(Ti)TiATj} 
DEM £ {DEM(Ti)ATj} ' 1 ' 
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Fig. 3. — Magnetic flux and loop temperatures for a sample of 
four bipoles. The black curves (left axes) correspond to the in- 
tegrated photospheric flux density of the bipoles. The red curves 
(right axes) are the emission weighted, and averaged coronal tem- 
perature profiles for the respective bipoles. Panels (a), (b), and (d) 
are for emerging bipoles, and panel (c) is for an evolved bipole. 

later for comparison with loop models. 

Integrated unsigned magnetic flux of both polarities as- 
sociated with these examples, as a function of time, are 
also calculated from HMfl Such profiles of temperature 
and magnetic field for a sample of four bipoles are plotted 
in Figure [3] Three cases of newly emerging bipoles (pan- 
els (a), (b), and (d)), and a case of an emerged bipole 
(panel (c)) are shown. The black curves are time profiles 
of magnetic flux (10 19 — 10 20 Mx), and red curves are 
temperatures in the range of 1 — 2 MK. 

Though they all fall in a category of emerging/ emerged 
loops, there is no clear relation between the magnetic 
flux at photosphere and the coronal loop temperature. 
In other words, it is not trivial to directly relate the field 
changes in photosphere to the temperature fluctuations 
in the corona. For example in Figure [31(a) , there is a 
strong correlation between the two physical quantities 
in the long term trend, but in panel (b) the tempera- 
ture seem to increase while the flux decreases. In pan- 
els (c)-(d), it is more complicated. We suggest that ev- 
ery emerging bipole may behave differently owing to its 
surrounding structures both in corona and photosphere. 
However, a common signature is that the temperature 
fluctuates/ rises at some stage in the emergence process. 

To further illustrate this behavior, we consider an- 
other example of an emerged bipole. In Figure |4] we 
plot the magnetic flux (black curve) and temperatures 
(red curves) this bipole. A sample image of this exam- 
ple from a particular time is also shown above (magnetic 
structure), and below (coronal loop) the plot. The thin 
vertical line demarcates the time of snapshot. The solid 
and dashed red curves are average temperatures derived 
from two adjacent regions (marked with solid and dashed 
lines in the image below the plot), of 2". 4 x 15" size each. 
In the accompanying animation, it is observed that, the 
drop in temperature after 500 minute is due to the re- 
connection (in the corona) of the parent bipole with the 



where AT!; = 0.1 in log (T), is the width of temperature 
bin around Tj. The time variations of Tdem will be used 



The two polarities are separated by a distance of approximately 
10 - 15 Mm in the photosphere. 
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Fig. 4. — Same as in Figure [3] Integrated flux density (black 
curve), and the temperatures (red solid and dashed curves) of an 
evolved bipole are plotted. A snapshot of the photospheric field 
configuration (top image, from HMI, 53" X 39"), and the corre- 
sponding coronal loop structure (bottom composite image, from 
AIA 171 A, 193 A, 211 A channels, 60" X 48") are shown for a par- 
ticular time as demarcated by the thin vertical line. The solid and 
dashed red curves are the emission weighted temperatures derived 
from 4 X 25 pixel (2". 4 X 15") rectangular boxes, from the regions 
marked in the bottom image respectively. A complete observed 
evolution of this example is presented as an animation, accompa- 
nying this figure. 

adjacent opposite polarity regions, changing the topol- 
ogy of the field, and completely disrupting the main loop. 
Hence the observed temperature of the loops originating 
from small ephemeral regions possibly depends on vari- 
ous factors. 

3. LOOP MODELING 

Temperature profile of the loop is a good diagnostic for 
the loop dynamics but to get a better picture, we also 
need to estimate the heating rate required to produce the 
observed temperatures. To this end, w e use Enthalpy- 
based Thermal Evolutio n of Loops (jKlimchuk et al.l 
[200l ICargill et al] [20121 ). EBTEL is a time depen- 
dent zero-dimensional (OD), hydrodynamic coronal loop 
model. For a given loop half-length and volumetric heat- 
ing rate, the code returns the loop properties in terms of 
average temperature, density, and pressure of the loop 
and also the values of these quantities at the loop apex 
(see the Appendix [AJ . 

We use EBTEL to model and derive the properties 
similar to the observed loops (we consider the example 
shown in Figure [31(a) for this purpose) . The properties 
include the DEM as a function of temperature, and the 
emission weighted temperature. We compare three dif- 
ferent heating scenario and discuss the results. For the 
models presented in next three subsections we make the 
following assumptions: (a) A loop is comprised of hun- 
dred individual strands, each with a constant length^ of 
about 18 Mm, and a uniform radius of about 0.1 Mm. In 
Figure [S]we plot the length and radius of a single strand 
as a function of time (thick solid and dashed lines, re- 
spectively), (b) Each strand is randomly heated with a 
certain heating profile over a period of 500 minutes, (c) 
The average values of various physical quantities over all 

3 This is only a rough estimate of the length based on the foot- 
point separation in the photosphere. 
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Fig. 5. — Physical dimensions of the strands used in this study. 
Thick solid and dashed lines, respectively, are the half length and 
radius of a constant strand. Thin solid and dashed curves, respec- 
tively, are the half length and radius of an expanding strand. 

the strands, represent the properties of the whole loop. 
Along with these assumptions, the heating events are 
chosen such that the modeled emission weighted temper- 
atures closely match the observed temperatures. 

In the Section 13.11 we describe the medium frequency 
heating model. Section [3.21 deals with the low frequency 
heating model. A medium frequency hybrid heating 
mod el is discussed in the Sections 13.31 and 13.41 In Sec- 
tion [33] we present an alternate explanation for the ob- 
served DEMs by considering a non-uniform cross-section 
of the loop. 

3.1. MEDIUM FREQ UENCY HE A TING MODEL 

In the medium frequency heating model (case 1), in- 
dividual strands are randomly heated with heating rates 
having 50 - 100 s temporal fluctuations. These rates are 
generated by a sequence of random numbers, and fur- 
ther filtering the signal within the desired band of peri- 
ods. The base or minimum heating rate is 10 -6 erg cm~ 3 
s , and the amplitude of the fluctuations vary by up to 
four orders of magnitude. The average heating rate for a 
single strand, over the entire duration of 500 minutes is 
about 4 x 10~ 3 erg cm~ 3 s _1 . 

A representative heating rate for one of the strands is 
plotted in Figure UKa). The plasma is reheated continu- 
ously before it is cooled to the equilibrium temperature 
due to base heating. In panel (b) we plot the resulting 
temperature of the strand apex (black, left axis) along 
with the resulting strand density (red, right axis). The 
temperature variations are 1 — 3 MK within in a single 
strand. Panel (c) is the average heating rate of all strands 
as a function of time. It should be noted that the fre- 
quency of this average quantity is not a relevant factor in 
distinguishing between various heating cases. Similarly 
in panel (d) we plot the average loop apex temperature 
(black, left axis), and the average loop density (red, right 
axis). Since the observed temperatures are derived from 
weighing the emission distribution, temperature of loop 
apex in panel (d) cannot be directly compared with its 
observed counterpart. 

In the top panel of Figure[7]we plot the observed DEMs 
(black dots). The observed DEMs have a broad distri- 
bution in temperature with a peak at log T of 6.3 and 
another peak at log T of 5.7 (same as Figure [3a)). The 
modeled DEMs, which have a narrow distribution, are 



5 




Fig. 6. — Results from a medium frequency heating model (case 
1). All the hundred strands are heated with an approximately same 
average heating rate, (a) A representative input heating given to 
a single strand with temporal fluctuations of 50 - 100 s. The base 
heating rate is 10 -6 erg cm -3 s — 1 for all the strands. The ampli- 
tudes of the heating rate fluctuate up to four orders of magnitude, 
(b) The resulting temperature of the strand apex (black, left axis), 
and loop density (red, right axis) for the heating profile shown 
in (a), (c) The average heating rate of hundred random realiza- 
tions, (d) The average loop apex temperature (black, left axis), 
and the average loop density (red, right axis) averaged over hun- 
dred strands, representing an observed loop. 

plotted as red dots with a similar temperature offset. 
Since there are heating events occurring almost contin- 
uously compared to the cooling time of the strands, the 
loop has no time to cool down completely and the tem- 
perature stays steady, with small fluctuations. Because 
of this reason, all the emission comes from a narrow dis- 
tribution of temperatures, which is reflected in the mod- 
eled results. 

Bottom panel is the resulting emission weighted tem- 
peratures from observations (black) and modeling (red). 
Note that the range, and level of fluctuations in the tem- 
perature match very well but, modeled DEM has a peak 
at log T of 6.25, and the predicted emission about this 
temperature is at least an order of magnitude more than 
the observed values. Further, the model predicts a weak 
or no emission at lower T. By increasing the magnitude 
of heating rate to match the temperature at which the 
peak emission occurs, will inherently increase the emis- 
sion, and also the weighted temperature well beyond the 
observed T. 

3.2. LOW FREQUENCY HEATING MODEL 

In the low frequency heating model (case 2), each 
strand is impulsively heated five times with an average 
of 100 minutes interval between each impulse. Each tri- 
angular pulse has a width of 500 s and a peak input of 
I0~ 2 erg cm~ 3 s" 1 . Further, the base heating remains 
the same as in case 1. In Figure Eta) we plot a sample 
profile of heat input given to one of the strands. Panel 
(b) is resulting temperature and density. Note that once 
the temperature reaches a maximum value, it takes about 
70 minutes for the strand to completely cool down. 

The average heating rate in panel (c) is less by a fac- 
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Fig. 7. — Comparison of DEM results obtained for case 1 with 
the observations. Top panel: DEMs from observations (black dots), 
and the modeled DEMs (red dots) are shown. All DEMs are given 
small temperature offsets for a better visualization of the distribu- 
tions. The black and red solid lines are the temporal medians of 
observed and modeled DEMs respectively. Bottom panel: Emis- 
sion weighted temperature derived from observations (black) and 
modeling (red). 




Fig. 8. — Same as Figure [6] Results from a low frequency (im- 
pulsive) heating model (case 2). On an average, each strand is 
heated every 6000 s once, with a triangular heating pulse having a 
maximum of 10 — 2 erg cm -3 s — , and a width of 500 s. 
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Fig. 9. — Same as Figure[7] but obtained for case 2. 

tor of 5 compared to that of case 1. It can be seen that 
the temperature (panel (b)) in this model has a broad 
distribution, which is reflected in a very broad DEM dis- 
tribution shown in Figure [5] (top panel, red). At log 
T of 6.5 the model produces a well constrained DEM 
that is higher than the observed DEM, although the ob- 
servations are less constrained at those temperatures. 
Also, the model predicts an overall higher emission at 
log (T) « 6.0 

The predicted emission weighted temperature (bottom 
panel, red) is comparable with observed temperature 
(bottom panel, black) . The level of fluctuations and the 
short term trend in the red curve are higher than what 
is seen in the observations. Furthermore, if the number 
of heating events are fewer than what is considered here 
(five), but with stronger impulses, the fluctuations now 
become noticeably large, and the observations should re- 
veal these features. 

3.3. MEDIUM FREQUENCY HYBRID HEATING 
MODEL 

For cases 1 and 2 we adjusted the model parame- 
ters such that the DEM- weighted temperature (Tdem) 
roughly matches the observed temperature for region 1. 
However, we find that the overall structure and features 
of the predicted DEM(T) does not match the DEM(T) 
derived from AIA observations. Therefore, neither of 
these models are fruitful in describing the 1 — 2 MK 
emerging loops in the quiet Sun. We suggest that the 
heating events may have a broad range, and/ or a popu- 
lation of different heating amplitudes, influencing differ- 




FlG. 10. — Same as Figure [6] Results from a medium frequency 
hybrid heating model (case 3). 20% of the strands are subjected to 
higher average heating inputs, but with same temporal fluctuations 
as in case 1 (see text for details). Shown in panel (a) is an example 
of lower heating rate case. 

ent strands. 

From the observational point of view, each strand in 
the loop is dynamically evolving and the lifetime of this 
unit is not clearly known. New strands emerge with the 
photospheric flux and replace the older ones in the loop. 
To investigate this problem further, a hybrid heating 
model (case 3) has been considered. In this model, we 
start with a simple assumption that 20% of the strands 
are rapidly heated with an excess amount of average 
heating rate of fifty times more than the remaining 80% 
of the loops. All loops receive a base heating similar to 
that of cases 1 and 2. 

In Figure [10] we plot the heating rates, temperature, 
and densities and also the respective average quantities. 
The profile in panel (a) is a low amplitude heating for 
a strand in the 80% population. The profiles shown in 
panel (b) are similar to that of Figure |B£b) , except for 
the overall lower values. Panels (c) and (d) show results 
averaged over all the strands, including the 20% that re- 
ceive a higher level of heating. In Figure [11] we show the 
DEM results for this case. The top panel is for the ob- 
served (black) and predicted (red) DEMs. We see that 
the predicted DEMs now have two distributions, clearly 
originating from the two populations of heating events. 
It is interesting to note how closely the observed and pre- 
dicted DEMs match. The emission weighted temperature 
is shown in the bottom panel of Figure 111! The fairly 
well reproduced quantities from this model suggest that 
a coronal loop, which has a bundle of many strands, can 
be heated by considering different amplitudes of medium 
frequency heating events. This is certainly a plausible 
assumption because, these emerging bipoles evolve con- 
tinuously, and various reasons can contribute to different 
heating episodes. 

Alternatively, we can also assume that each strand 
spends 20% of its time being heated to higher values 
(similar to the 20% strands case described previously), 
and the remaining time to lower values. Both the see- 
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Fig. 11. — Same as Figure [7] but obtained for hybrid heating 
model (case 3). 




Fig. 12. — Alternative version of case 3 in which, each strand 
receives high heating events 20% of its time, and the remaining 
time, low heating events. In panel (a), a sample heating profile is 
plotted. The panels (b) shows strand apex temperature (black, left 
axis) and strand de nsit y (red, right). The panel (c) is same as the 
top panel of Figure [TT] 
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Fig. 13. — Hybrid heating model with expanding length and ra- 
dius as shown in Figure \5\ (thin solid and dashed curves respec- 
tively). Both the length and radius of the loop vary slowly with 
time. In this case the temperature (red, lower panel) increases as 
the length of the loop increases. 

narios produce similar results. In Figure [12] we plot the 
results from this alternative case. The panel (a) is a 
sample heating profile showing both low and high heat- 
ing events. In panel (b) the loop apex temperature and 
density are shown. In panel (c) the observed (black), 
and predicted (red) DEMs are plotted along with their 
respective temporal medians. 

3.4. A CASE OF EXPANDING LOOP 

In general, the coronal loop length increases with time 
as it emerges through the solar atmosphere. Also, the 
area as a whole, as the strength of the magnetic field 
drops with height, the area of strand increases with time. 
Due to this expansion, filled in plasma may experience 
additional adiabatic cooling effects, as the loop pres- 
sure and density are modified by the volume change (see 
the Appendix [Aj . We tested heating model described in 
case 3 on a slowly expanding loop, comprised of hundred 
strands as explained in the previous sections. 

The half length and radius of a single strand are shown 
as thin solid, and dashed curves in Figure [5] The DEMs 
and T are plotted in Figure [13] The way we consider the 
volume expansion is that, each strand slowly expands for 
about 350 minutes and then the expansion saturates to 
a constant value. This constant value matches with the 
length, and radius of the strand chosen in all the cases. 
It is observed that for a slowly expanding loop, the adi- 
abatic cooling effect can be negligible. There are two 



competing effects here. Under equilibrium conditions, 
the temperature of the loop increases with the length. In 
our slowly expanding loop, the cooling is compensated for 
with the length increase. But in reality, the rate of vol- 
ume expansion can be entirely different, and more rapid 
than what we considered here. These effects become im- 
portant when changes in the loop pressure and density 
due to expansion alone, and heating are comparable. 

3.5. EFFECTS OF NON- UNIFORM 
CROSS-SECTION OF THE LOOP 

In the above EBTEL-based models the emission is as- 
sumed to come from the coronal portion of the loop. 
EBTEL also predicts the emission from the transition 
region (TR) to model the lower temperatures. However, 
the predicted TR emission is strong and rather flat rela- 
tive to the corona. Inclusion of up to 5 — 10% of the TR 
emission will not affect the results, but adding more con- 
tribution from the TR requires stronger heating to match 
the observed emission weighted temperature. This re- 
sults in a strong emission from the higher temperature, 
which is not observed. 

One possible reason for the strong TR DEMs produced 
by EBTEL is that the model assumes a constant cross- 
section over the length of the loop, whereas the loops 
on the Sun have significant expansion factors (7) be- 
tween the loop footpoints in the TR and the loop top 
in the coro na. Potential field modeling of active re- 
gions (e.g., lAsgari-Targhi fc van Ballcgooiien 2012) in- 
dicates 7 = 3 — 30, depending on height, and similar ex- 
pansion factors may occur on the quiet Sun. When the 
cross-sectional area A of a loop increases with height, the 
volume of plasma at coronal temperatures is increased 
relative to that at TR temperatures, so the slope of 
the DEM(T) curve becomes steeper and more consistent 
with observations. 

To demonstrate this effect, we developed a simple loop 
model for the case that the cross-section A varies along 
the loop. The heating is assumed to be steady in time. 
The model is described in Appendix [B] It allows us to 
compute the DEM(T) for a single loop with a given ex- 
pansion factor 7, half-length L, and peak temperature 
T max (we use L = 9 Mm). We repeat the calculation 
for different peak temperatures (6.0 < log 10 T max < 6.4) 
and compute the average DEM(T). In Figure Q3] we 
plot the DEM results for 7 = 1 (left panel) and 7 = 5 
(right panel). Note that for a loop with uniform cross- 
section (7 = 1), DEM(T) is flat for T > 10 5 K, simi- 
lar to the DEMs pre dicted with the EBTEL code (see 
iKlimchuk et al.l 120081 ). In contrast, for 7 = 5 the peak 
value of the DEM in the corona is about 8 times its value 
in the TR, similar to the observed DEMs (see Figure [2]). 
These results suggest that the overall shape of the ob- 
served DEM can be very well reproduced with a collec- 
tion of hot loops (T max > 1 MK) that have significant 
expansion factors (7^5 -10). However, the peak value 
of the DEM as predicted by the model is larger than the 
observed value by a factor of about 100. Therefore, the 
loops must fill only a small fraction of the coronal volume 
(filling factor ~ 1 %). 

4. SUMMARY AND DISCUSSION 

Using the high temporal cadence observations from the 
HMI and AIA instruments on board SDO, we studied 



the cases of emerging bipolar regions in the quiet Sun. 
High cadence data from AIA including six EUV channels 
are re-sampled to 1 minute data to improve the signal-to- 
noise ratio, as well as, to have a good temporal resolution. 
Further, xrt_dem_iterative2 .pro is used to construct 
DEMs near loop top in a 6" x 6" pixel region (Section El 
Figured]). From these DEMs, we get the temporal evolu- 
tion of emission weighted temperature with Equation ([2]) . 

Integrated unsigned magnetic flux derived from the 
HMI observations (10 19 — 10 20 Mx) is compared with 
the temperature of the loop, for each example (Figure[3|). 
There is no clear relation between the two quantities, sug- 
gesting that, for these small emerging bipoles, the sur- 
rounding regions in photosphere and higher atmosphere 
play an important role in the loop evolution. 

To estimate the energetics involved in the formation of 
these loops we use a hydrodynamic loop model (EBTEL) 
to simulate the DEMs and emission weighted tempera- 
tures. We assume that a loop is a bundle of one hundred 
strands, each having a length of about 18 Mm, and a 
uniform radius of 0.1 Mm. Furthermore, each strand 
is randomly heated and the average effect describes the 
properties of the observed loop. To this end, we tested 
three simple heating events with varied heating frequen- 
cies as described in Sections 13.11 13.21 13.31 and 13.41 The 
average heating input in our study (» 10 6 ergs cm -2 s _1 ) 
is in close agreement with the approximate energy losses 
observed in the quiet Sun. The 3cr values of the fluctua- 
tions in To em (MK) are about 0.25, 0.16, 0.55, and 0.24 
for the observations, cases 1, 2, and 3(a) respectively. 

In case 3, we tested the following sub-cases: (a) 20% 
of the strands are heated to high heating values all the 
time, and the remaining strands are heated to low heat- 
ing values (Section l3.3|) . (b) all strands are heated to high 
heating values 20% of their time (Section l3.3[) . (c) similar 
to case 3(a), but for expanding str ands to account for the 
adiabatic cooling effects (Section |3.4[) . It is shown that 
the cases 3(a), and 3(b) are equivalent and match the ob- 
servations fairly well. This suggests that there may be a 
range of heating events operating in the loops at a given 
time. In the cases 2 and 3(b), though there are only a 
few large heating events, the essential difference between 
the two cases is that, unlike in case 2, the duration of a 
single high heating phase in case 3(b) itself is longer com- 
pared to typical plasma cooling time (making case 3(b) 
statistically a steady heating model). This allows case 
3(b) to find a DEM peak at higher temperatures. The 
model described in case 3 is the best model we could 
obtained with in the scope of the present work. Mix- 
ing low and medium heating at various proportions with 
different average heating rates show discrepancy, and do 
not fit observations completely. These results are based 
on the assumption that the emitting plasma has coronal 
origin. 

Alternatively, we also argued that to include TR emis- 
sion in the model, it is important to consider an expan- 
sion of the loop from TR to corona. In this scenario, a 
steady heating model for loops with loop apex temper- 
ature > 10 6 K can well reproduce the observed DEMs, 
assuming ~ 1% plasma filling factor. 

Rel iability of AIA DE Ms is a matter of de- 
bate. lO'Dwver et al.l (|2010T ) studied the contribution of 
spectral lines and continuum emission to the AIA EUV 
channels using CHIANTI atomic database. They em- 
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Fig. 14. — Results from a simple loop model assuming steady heating. Left panel: DEM of a loop with expansion factor (7) = 1. Note 
the fiat DEM in the rang e of 10 5 < T(K)< 10 6 . Right panel: Similar to the left panel, but with 7 = 5. The DEM now shows a clear peak 
atTR3 2 MK. 



phasize that the contribution of particular spectral lines 
and continuum emission can affect the interpretation of 
the observed features, when AIA channels are used to 
observe regions other than those for wh ich the channels 
were designed. IDel Zanna et al.l (|2011A compared AIA 
DEMs with the HinodefElS observations of active re- 
gions. They found discrepancies between the derived 
DEMs. This is mainly due to the multi-thermal na- 
ture of AIA response curves, which have contributions 
from cooler components. The cooler emission below 6.0 
(in log T), seen in our observations could be due to 
th e double-peaked nature o f AIA responses as suggested 
by IDel Zanna et "ail (|201 If ) . Empirically modified filter 
response curves for AlA are derived to address some of 
these issues however, the role of this possible contamina- 
tion in a already cool loop (like the one originating from 
a small bipole in the quiet Sun), as compared to the 
warm loops in the hotter active region has to be further 
examined. 

The models presented in this work assume that the 
strands are heated uniformly over their entire length. Al- 
ternatively, the strand can be heated in a non-uniform 
manner with localized and concentrated heat sources. 
If the heating is concentrated at the loop footpoints, 
this may lead to the loss of equilibrium in the energy 
balancing terms, as the radiative losses in the coronal 
section dominate the downward conductive flux. This 
will trigger the runaway cooling due to strong radia- 
tive losses and a con densation is formed in the coro- 
nal lo ops (for example [Hood fc P riest 1980t iMiiller et all 
2004). This is a well studi ed phenomenon in the forma- 
tion of solar prominen ces ()Antiochos &: Klimchuklll991l ; 
lAntiochos et al.l 119991 ) . Recently, based upon the ob- 
served pr operties of the hot, an d warm loops in active 
regions, Klimchu k et al.1 ((2010) have argued that the 
high concentration of heating low in the corona, and 
the steady or quasi-steady heating models (leading to 
thermal nonequil ibrium) can be ruled out. However, 
I Peter et al.l (|2012l ) claim that a steady supply of energy is 
required even in the events of condensation in the corona, 
to keep the coronal pressure. They also suggested that 
thermal non-equilibrium can be a valuable tool in in- 
vestigating the plasma dynamics and heat input in the 



regions where condensation forms. 

The studies on the role and importance of the thermal 
non-equilibrium in the formation of conde nsation in the 
short q uiet Sun loop s are not extensive. IMuller et al.l 
l|2003f) . IMuller et ail (|2004D discussed in detail the nu- 
merical simulations of condensation and catastrophic 
cooling of short TR 10 Mm loops, and longer 100 Mm 
coronal loops, respectively. They considered heating that 
has exponential height dependence along the loop, and 
further suggested that the catastrophic cooling is initi- 
ated by the loss of equilibrium at the loop apex due to 
concentration of heating at the footpoints, but not due 
to a drastic decrease of the total loop heating. 

Note that the strands in a loop may interact in a very 
complex manner, and their response to the condensation 
is the key objective to be addressed. Further work is nec- 
essary to get a better picture of the nature and location 
of the heating, observational signatures of condensation, 
and finally the role of magnetic field in this whole pro- 
cess. A complete set of answers for these questions is 
still elusive and we need more observational constraints 
to narrow down the possibilities. 
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APPENDIX 

EBTEL AND LOOPS THAT EXPAND WITH TIME 

The standard version of EBTEL assumes a symmetric loop with constant loop length and, uniform cross-section. 
The model is based on the ID time-dependent energy conservation equation 

f = -^ + p)-^ + o-» 2 A (n (ad 

where s is a spatial coordinate along the magnetic field; E = |P + \pv 2 is the total energy density; n, T, P, and v are 
the electron number density, temperature, total pressure, and plasma bulk velocity, respectively; F c is the heat flux; 
Q is volumetric heating rate; and A(T) is the radiative loss function for optically thin plasma. It is assumed that the 
velocity and heat flux both vanish at the loop apex due to symmetry. Also, the flow velocity is subsonic, and gravity 
is neglected in the energy equation. Integrating the above equation over the coronal (L), and TR (I) lengths with the 
above assumption, we get 

3 dP 5 

-L—n-P v +F +LQ-Tl c , (A2) 
3 BP 5 

~l-^-™--P vo-F + lQ tr -nt r , (A3) 

where overbar denotes the spatial averages of the quantities over respective sections of the loop, and subscript 
denotes the values at the base of the corona; 7Z C , and lZt r are the coronal, and the TR radiative loss rates respectively. 
Neglecting the terms involving I (for a thin TR), and together with ideal gas law0, P and n can be approximated^ 
with 



dP ~ 2 
~dt ~ 3 

dn 1 



Q-^{n c + n tr 



dt 5kLTo 



(A4) 

(Fo + Ktr). (A5) 



For a given heating rate Q(t), the EBTEL model returns P, n, and T, with other useful quantities. 

For a uniformly expanding strand of length L(t) and radius R(t) adiabatically, the above equations are modified by 
adding a term — jP£(t) on the right hand side of Equation (|A4I) . and — n£(t) in Equation (|A5p . where 7 = 5/3 is the 
ratio of specific heats, and 

1 dL 2 dR ,,„. 

= (A6) 

The strand pressure, density, and temperature are modeled accordingly. The time varying length and radius of a single 
strand are shown as thi n sol id and dashed curves, respectively, in Figure [S] Note that R(t) explicitly enters the scheme 
only through Equation (IA6|) . and everywhere else, it is absorbed due to volumetric averaging. 

MODEL FOR LOOPS THAT EXPAND WITH HEIGHT 

In this section we describe a loop model for the case that the cross-sectional area A varies along the loop. For 
simplicity the area A(T) is considered to be a function of temperature: 



A(T) = exp\ In 7 



(Bl) 



ATn: 

where T max is the maximum temperature at loop top, and z(T) is a monotonically increasing function, starting with 
z « at the base of the TR. We use z(T) = y + \J\ + y 2 with y = (x — xq)/x\ and x = log 10 T, The constants £0 and 
X\ are set to 6.0 and 0.2, respectively, so that most of the area change occurs near a temperature of 1 MK. Similarly, 
the volumetric heating rate is 

Q(T) = Q max (7^—) , (B2) 

\ ^ max / 

where Qmax is the heating rate at the loop top, and m is an exponent (for the models presented here we set m = 0). 
The loop is assumed to be symmetric, and heating is assumed to be steady in time. We solve the following energy 
balance equation: 

^ (AF C ) = A(T) [Q(T) - n 2 A(T)] , (B3) 

4 P = 2nkT, where k is the Boltzmann's constant. 5 ■ , , . ■ . . , , . ., 

simple volumetric averaging yields similar results. 
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where s is a spatial coordinate along the magnetic field, F c (s) = —K^T b / 2 dT / d s is the conductive hea t flux, n(s) is 
the electr on d ensity, and A(T) is the radiative loss function, which is taken from iKlimchuk et al.l (|2008| ). Multiplying 
equation (IB3|) by AF C and integrating over position along the loop, we obtain: 

\A 2 F 2 = K Q ma MEi(T) - E 2 (T)}, (B4) 

where 

Ex(T)= [ A 2 (T)A(T)T 1/2 dT, (B5) 

J Tbaso 

£ 2 (T) = / A 2 (T)(T/T m ^) m T 5 / 2 dT. (B6) 

Here Tb aS e is the temperature at the base of the TR (Tbasc = 10 4 K), and we assume F c = at the base. The factor / 
is given by 

f = p2 = ^ 2 (^" lax ) (-R7\ 
~~ 4k 2 EifT V 

where P = 2nkT is the plasma pressure (a constant), and the last equality in (|B7|) follows from the requirement that 
= at the loop top. Then the loop half-length L is given by 

For a given peak temperature T max and half-length L, we can compute the heating rate Q max , pressure P, heat flux 
F C (T), and density n(T). Then the DEM is given by ip(T) = n 2 {T)A{T){dT / ds)- 1 . Since the area factor is normalized 
such that A(T max ) = 1, this DEM(T) does not include the effects of a possible filling factor of the coronal loops. 
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